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A unified (classical-quantum-statistical) formalism for a system with continuous spectrum is 
introduced. For this kind of systems ergodicity behavior and the existence of microcanonical and 
canonical (KMS) equilibrium is proved. It is argued that the continuous spectrum condition is 
essential for the thermodynamical behavior. 
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O ■ I. INTRODUCTION. 

O ; 

I Systems with continuous spectrum have been studied at the classical, quantum, and statistical level. In papers [Q, 
0, and || we have contributed to this research. Precisely, based on the ideas of van Hove Q we have introduced a 
method to deal with these systems, that we have successfully used to treat decaying phenomena, statistical mechanics 
problems and to find the classical limit of quantum systems ||. In this paper, using the same technique, we would 
like to present a panoramic version of the different "limits"]] between physical levels. We will see that we can use our 
00 ' method (complemented by other ideas) as a unified formalism that allows to go from one level to the other. Even if 
some features of the problem are well known from a long time we believe that a unified treatment of these matters, 
like ours, was missing. Furthermore, most probably, it will be quite useful. Essentially we believed it is useful to see 
how the interplay of classical and quantum concepts can be used to solve problems. 
So let us begin looking at the whole panorama. Physics studies two kind of systems: 

A. - Individual systems, like particles or single physical systems. This is the case in classical mechanics, electromag- 
netism, special, and general relativity. 

B . - Statistical systems, like ensembles or sets with many copies of a single physical system. This is the case in 
quantum mechanics, statistical (classical or quantum) mechanics, thermodynamics, and quantum field theory. 

In category A the limits between different levels are simple and well known: 
Ai.— We can go from general relativity to special relativity if we consider just flat space-times. 
A2.— In the "classical limit" (/? — - — > 0) special relativity becomes classical mechanics. 
The limit between the two categories, e. g.: of statistical mechanics towards classical mechanics is produced by 
the phenomenon of localization that we have study in paper ||, where the necessary conditions under which this 
phenomenon take place (localizing potential, adequate initial conditions, etc.) are given for some examples^. So we 
will just make a short comment of this problem in the paper. 
q-( In category B the limits are much more involved and will be considered in this paper. So it is organized as follows: 
Bi.— The thermalization of a quantum system is considered in section II. We will study it in two steps: 
Bi.x-— A first step is decoherence (section II. A). This section, as section III, it is a brief resume of ref. |J included 
in the paper for the sake of completeness. The conditions under which the evolution of the quantum system decohere 
5_] ' (more than one bound state, small decoherence time, etc.) are listed in the just quoted paper 
B1.2.— We will see how we can obtain a KMS equilibrium after decoherence (section II. B). 
B2.— The quantum to classical statistical mechanic limit is studied in section III. . 

B3 _ The statistical classical mechanics relation with the theory of dynamical system and with thermodynamics 
will be considered in two steps in section IV: 

B3.1.— In the first step the classical statistical system is considered as dynamical one. Then the micro-canonical 
ensembles appear as a consequence of dynamical systems theorems (section IV. A). 

B3.2.- Finally, in the second step, the canonical ensembles and thermodynamic are introduced (section IV. B). 

Other limit can be obtained combining those above. Knowing all these matters we can foresee the final fate of a 
quantum system (to remain a quantum one, or a classical one, or to end in thermodynamical equilibrium, etc.). 



1 Usually the relations between physical levels are called "limits", even though only in some cases they are just simple mathe- 
matical limits. 

2 The study of the localization phenomena is a deep and complicated subject, since in its classical version is related with the 
problem of the "limit circle" (see |^|, cap II). This is one of the problems that Hilbert listed for the XX century mathematicians 
and it remains still unsolved. 
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Section V is devoted to the localization problem (namely the B — > A limit). 
We will state our conclusion in section VI. 



II. QUANTUM TO STATISTICAL LIMIT. 



We will demonstrate that for a wide set of quantum systems non-diagonal density operators can be considered as 
the transient phase while diagonal density operators is the final regime and a permanent stateQ. We will find a basis 
where exact density operators decoherence appears for these final states . Then we will see how the quantum systems 
thermalize. 



A. Decoherence. 



1. Decoherence in the energy. 



Let us consider a closed and isolated quantum system with N + 1 dynamical variables and a Hamiltonian endowed 
with a continuous spectrum and just one bounded ground state. So the discrete part of the spectrum of H has only 
one value lo$ < ^ and the continuous spectrum is, let say, < d < oo. Let us assume that it is possible to diagonalize 
the Hamiltonian H, together with N observables Oi (i = l,...,N). The operators (H, Oi,...,On) form a complete 
set of commuting observables (CSCO). For simplicity we also assume a discrete spectrum for the N observables Oi. 
Therefore we write 

POO 

H = ujq J]] 1^07 m ){ ljJ o^ m \ + / w E] |w, m) (ui, m\du (1) 

m m 

where l>o < is the energy of the ground state, and m = {mi, ...,itin} labels a set of discrete indexes which are the 
eigenvalues of the observables 0\,...,On- {|woi m )i |w,rn}} is a basis of simultaneous generalized eigenvectors of the 
CSCO: 

H\luq, m) = wo|wo, m), H\u),m) — u)\u>,m), 
Oi\ujQ,m) = mi\u>o,m), Oi\u,m) = mi\u,m). 

As at the statistical quantum level we only measure mean values of observables let us define the observables we will 
use. The most general observable that we are going to consider in our model reads: 

/•oo 

0=y2 O(u )mm'\^o,m){u ,m'\ + V] / duiO(ui) mm >\u),m){oj,m'\ + 
i i Jo 

mm' 

dujO{uj,ujo) mm '\u,m){uJo,m!\ + 
du>'O(u ,w') mm >\uJa,m)(u;',m'\ + 

duidcj 0(u),uj ) mm i\LJ,m}(u) ,m |, (2) 



■E 

mm' 

E 

mm' 

r>oo poc 



E 

mm' 



where = O and O(oj) mm i, O(u,uo)mm'i O(wo,w) mm ' and 0(u),u)') mm i are ordinary functions of the real variables 
li and d' (these functions must have some mathematical properties in order to develop the theory; these properties 



3 The arguments of this section are fully developed in paper Q, where the philosophy of the method is explained at large. 
4 The case of many bound states is considered in paper Here we consider the case with only one bound state because in 
this case the system usually decoheres. 
5 Other quantum measurements imply a limit from category B to category A (cf. 0|) which is outside the scope of this paper. 
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are listed in paper Q). We will say that these observables belong to a space O (which is contained in the algebra 
studied in 0]). This space has the basis {\luq, mm'), \uj,mm'), \oj ujq, mm'), \loqlo' ,mm'), \uiui' ,mm')}: 



\u)o,mm') = \ujq, m) (loq, m'\, \uo,mm') = \uj,m)(uj,m'\, 



\ujLUo,mm') = \uj,m)(uJo,m'\, \loqlo' ,mm') = \luq, m)(uj', m'\, 



(3) 



\luuj' , mm') = \lu, m)(uj', m'\ 



The quantum states p are measured by the observables just defined, computing the mean values of these observable 
in the quantum states, i. e. in the usual notation: (0) p = Tr(p^O). We can consider that mean values are the more 
primitive objects of quantum theory 0. These mean values, generalized as in paper ||, can be considered as linear 

functionals p (mapping the vectors on the real numbers), that we can call (p\0) (see also [8]). Then p G S C O , 
where S is a convenient convex set contained in O , the space of linear functionals over O [9], |n]|]. The basis of 
O' (that can also be considered as the co-basis of O) is {{ujQ,mm'\, (uj,mm'\, (luloq , mm' \ , (luquj' ,mm'\, (ww',mm'|} 
defined as functionals by the equations: 



(uj ,mm'\uj ,nn') = 6 mn S n 



(lu, mm'\r\, nri) = 5{u - T])S mn dr, 



(ujujQ,mm'\riu! ,nn') = 5(u> - r])5 mn S mlnl , 



{u u',mm'\uj r]',nn') = 5(u>' - r]')5 mn S m > n/ , 



(lulu', mm!\r\rj ', nn) = 6(u> — rj)6(u>' — rj')8 mn 5r, 
and all other (.|.) are zero. Then, a generic quantum state reads: 



(4) 



p OG 

P=T] K^o) mm ,(u>o,mm'| + V / du)p(u) mm/ (u,mm\ 
, i Jo 

mm' mm' 
roc 

+ V / du;p(uj,u;o) rn7n/ (Luu;o,mm f \ + 

mm' 

+ V / daj'p{u ,Lj') mm ,(uJouj',mm'\ + 
, Jo 



E 



du I du' p(uj,uj') mm ,(ujLu' ,mm'\ 
o Jo 



(5) 



where 



p(LO,UJ ) mm , = p(LU( ) ,Uj) m > rn , p{LO,u') mm , = p{w',U)) m r m , 



and p(LOo) mm and p(u>) are real and non negative satisfying the total probability condition 

/>OC 

(p|J) = ^ p(vo)mm + E / duip(uj) mm = 1, 



(G) 



where / = ^2 m \loo, m) [luq, m\ + J Q °° duj ^2 m \lo, m) {lo, m\ is the identity operator in O. Eq. (^) is the extension to state 
functionals of the usual condition Trp^ = 1, used when p is a density operator. Thus, from now on, Trp = (p\I)- 
The time evolution of the quantum state p reads: 



3 We can choose more general O but the one we have defined is enough for our purpose. 
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p oo 

p(t) = V pW^/h^m'l + V / rfwp(cj) mm/ (a; ) mm / | + 
/ / jo 

mm' mm' 
poo 

+ / ^K^) w e i ^--°) i (a; W o,mrn / | + 
/ Jo 

mm 

/•OO 

+ £ / ^'^M) mm -e i ( u »- w ') t (^>m'| + 

mm' 

P OO oo 

+ / <W ^'pM^^^'l'H'.mm'l (7) 

mm' 

As we have already said at the statistical quantum level we can only measure mean values of observables in quantum 
states, i. e.: 

(0) p(t) = (p(t)\0) = 

poo 

= V p(uo) mrn/ 0(cJ )mm> + V / dup(u)) mrn ,0(u)) m m' + 

/ / JO 

mm mm 

poo 

+ V / dL>p(u>, w ) mm ,e« u - u °) t O{u, wo)w + 

mm 

poo 

+ V / d W '^y) mm ,6 i ( U °- W '' t O(Wo ) W , ) mro - + 

/•OO pOO 

+ ]T / dw / dw%^ mm reM t O(«,w')m«,' ) (8) 



.,.., ,Jo Jo 

Let us now consider the fate of these mean values when t — ► oo, using the Riemann-Lebesgue theorem we obtain the 
limit, for all O € O 

lim <0)„ w = (0} p , (9) 

t— >oo V ' 

where we have introduced the diagonal asymptotic or equilibrium state functional 

p OO 

P* = P( w o)mm'( w o,mm'\ + V / dujp(uj) mm ,(uj,mm'\ (10) 
/ / Jo 

mm' mm' 

Therefore, in a weak sense we have: 

W lim p(t) = p* (11) 

t — >oo 

Thus, any quantum state weakly goes to a linear combination of the energy diagonal states (ijjQ,mm'\ and (u,mm'\ 
(the energy "off-diagonal" states {u)ujQ,mm'\, (u>qUj' , mm' \ and (ww',mm'| are not present in p*). This is the case if 
we observe and measure the system evolution with any possible observable of space O (albeit the discussion in section 
IV. A. 3). Then, from the observational point of view, we have decoherence of the energy levels, even that, from the 
strong limit point of view the off-diagonal terms never vanish, they just oscillate, since we cannot directly use the 
Riemann-Lebesgue theorem in the operator equation (^). 

2. Decoherence in the other "momentum" dynamical variables. 

Having established the decoherence in the energy levels we must consider the decoherence in the other dynamical 
variables Oi, of the CSCO where we are working. We will call these variables "momentum variables". As the 



7 The three of diagonal terms of eqs. (Q) or (g) disappear. This happens because we have just one bound state and would not 
be the case for more than one bound state. 
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expression of p* given in eq. ( |io| ) involve only the time independent components of p(t), it is impossible that a 
different decoherence process take place to eliminate the off-diagonal terms in the remaining N dynamical variables. 
Therefore, the only thing to do is to find if there is a basis where the off-diagonal components of p(too) mm ' and 
p(u)mm' vanish at any time. 

Let us consider the following change of basis 

\ UJ o,i') = ^C/(w )mr|w ,m), \oj,r) = y^f7(cj) mr \w,m), (12) 



where r and m are short notations for r = {r\, r^} and m = {mi, m^}, and [f/(f2) ] = U(fl) rm (fl denotes 
either Wq < or u> € R + ). 

The new basis {|o;oj r )j l w i r )} verifies the generalized orthogonality conditions 

(ujQ,r\ujQ,r') = S rr >, (w,r\u}',r ) = S(u> — ui')5 rr >, 
(u)q, r\ui, r') — (uj, r\u>o,r') = 0. 



As p(u>o) mm > = p(wo)m'm and p{to) mm , — p(ui) m i m , it is possible to choose U(u>o) and U(ui) in such a way that the 
off-diagonal parts of p(Lu ) rr i and p(uj) rr i vanish, i.e. 

p(tjJo)rr> = Pr(^o) $rr' , p{w)rr' = Pr{u) 0~rr> • (13) 

Therefore, there is a final pointer basis for the observables given by {\u>o,rr'), \u>, rr'), \ujujQ,rr'), \ujqlo' , rr' ) , \lulo' ,rr')} 
and defined as in eq. (|^). The corresponding final pointer basis for the states {(too, ( w i rr '\i (^luq, rr'\, (ujqlo' , rr'\, 
(ujlo' , rr'\{ diagonalizes the time independent part of pit) and therefore it diagonalizes the final state p* 



t — >oo 



oo 



p* = W )im p(t) = p r (oJo)(LJo,rr\ + ^ / dojp r (uj){uj,rr\. (14) 



Now we can define the final exact pointer observables ]l2|]P| 

/>oo 

Pi= VP;M|w ,r)(a;o,r|+ / dwV;j?(a;)|w,r>(w,r|. (16) 

As 7? and Pi are diagonal in the basis {l^o, r ), |w, r)}, the set Pi, ...Pat} is precisely the complete set of commuting 
observables (CSCO) ^ related to this basis, where is diagonal in the corresponding co-basis for the states. For 
simplicity we define the operators Pi such that P^(u)o) = P?(v) = r i, thus 

Pj|wo,r) = r t \u) ,r), Pi\w,r) = n\cj,r). (17) 

Therefore {\(-Oq, r), \oj, r)} is the observers' final pointer basis were there is a perfect decoherence in the corresponding 
state co-basis. Moreover the generalized states (loq, rr\ and (oj, rr\ are constants of the motion, and therefore these ex- 
act pointer observables have a constant statistical entropy and will be " at the top of the list" of Zurek's " predictability 
sieve" @. 
Therefore: 

L- Decoherence in the energy is produced by the time evolution. 



3 If we would like to just have integrals and treat all spectra as continuous we would write 

p oo p oo p oo 

^2 duj = / ... / dujdR 1 ...dR N S(r 1 - R 1 )...S(r N - R A 



j du j dp(R) (15) 
9 In usual quantum mechanical system this set is numerable. 
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ii.- Decoherence in the other dynamical variables can be seen if we choose an adequate basis, namely the hnal 
pointer basis. 

Our main result is eq. (pit): When t — > oo then p(t) — > p* and in this state the dynamical variables H,P\, ...,Pn 
are well defined. Therefore the eventual conjugated variables to these momentum variables (namely: configuration 
variables, if they exist) are completely undefined. 

In fact, calling by Lj the generator of the displacements along the eventual configuration variable conjugated to Pj, 
we have (hip*\0) = (p*\h\0) = (p*|[Pj, O]) = for all O e O as it can be proved by direct computation 0. Then p* 
is homogeneous in these configuration variables. 

From the preceding section we may have the feeling that the process of decoherence must be found in all the physical 
systems, and therefore, all of them eventually would become classical when H — ► 0. It is not so, this evolution only 
happens under certain conditions, as explained in ||. I. e., if there is more than one bound state or the decoherence 
time is infinite the system does not decohere. 



B. Quantum thermalization. 



1. The minimally biased state. 

At this point we can find the final quantum minimally biased state according to information theory Jl3| ] . For 
simplicity let us just use one observable in our CSCO, the hamiltonian H with no bound state Wo0- Then eq. ( |l4| ) 
reads 



{p,\ = I p(Lu){L0\du (18) 

(p*| has trace one, so from eq. (0) we know that 

Trp* = (p*\I)= I p{oj)doj = l (19) 



The mean value of the energy reads 

/•OO 

(H) Pt = (p*\H) = / up{u)dLo = E (20) 
Jo 

If this is the only data available to obtain the minimal biased state we must maximizes the (Shannon) missing 
information P| 



f-OO 

I = H[p*} = -K I p(uo) log p{u>)du) 
Jo 



(21) 



contrained by eqs. (19) and (poj) . The solution of this variational problem is p(cu) ~ e ^ so the unbiased is 

poo 

(p»| = Z- 1 / e-^(u\duj (22) 



where the constants Z and /? can be computed from eqs. ([l9]) and (|2p|). 

The next problem is to represent this minimal biased (p*| as e _/3 . But an equation like (p^l—e^ 1311 is impossible 
since H is a "vector" (ket) and (p*| is a "functional" (bra). To solve this problem we must go back to eq. ( [l9| ) and 
see that the trace of the state can be written as 

Trp* = (p|/) = (p| / \u)du; \I) = / \u)du (23) 



10 In fact, in our model H is the only member of the CSCO with continuous spectrum, so all the problem is contained in H. 
I. e., the continuous version of eq. (2.3) or Q or eq. (5.6) of jwj, postulating that when no information or probabilities 
are available we should consider the cells equally probable. 
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Symmetrically we can define the trace of an operator A as 



TrA=(I\A)= (Lj\du\A); (I\ = / (cj\duJ (24) 
Jo Jo 

We will see that this definition is completely reasonable at the end of section IV. B. Using it we can write the functional 
that corresponds to (p» | as 

W0 [A] = {pAA)= K -^ m j ( 25 ) 

We can immediately see that the trace of this functional is one since wp\I\ = {p*\I) = 1- Moreover, the last definition 
coincides with the one of eq. ( p2] ) since e 

-0H = j e-Pu^du and 

Z = {I\e~ 0H ) = J J {w\e- pu \u')dudu' = J e- 0w 8(u - uj')dujduj' = J e~ pu du (26) 

Also 

(I\e- 0H A) = f e-P"(Lj\A)du (27) 



From the two last equations it is evident that (|22|) coincides with fl25| ) and we can write the minimally biased canonical 
final state as 

{pA=Z-\l\e-P H (28) 
It is easy to generalize the above reasoning from H to the CSCO {H, Pi, ...P/v}. Then we would obtain 

(jh\ = z- 1 (I\e- f,H -i lP i---* irPtt (29) 
that corresponds to a generalized grand-canonical ensemble. 

2. KMS equilibrium. 

As a demonstration that t0/3|yl] is the good equilibrium limit at temperature /3 _1 we will prove that it satisfies the 
KMS condition Jl6] ]. The time evolution of operator A is given by 

a t (A) = e lHt Ae- lHt (30) 

( eq. (Q), for the state evolution, is deduced from this equation). Then we must check the analyticity properties of 

(I\e-P H a t (A)B) (I\ e -^e^Ae~^B) 
MMA)B] = {ile _ pH) = ^— ^ (31) 

We must first prove that definition (Q) has the cyclic trace property. Since the algebra O is associative ||] we only 
need to demonstrate the commutativity. In factf^ 

Tr(AB) = (I\AB) = ( J (w\du\AB) = J A^B^du = (I\BA) = Tr(BA) (32) 



then 



W/3[a t {A)B\ = (/| e -/3g) = (I\e-P H ) = w P{ Ba t+iP\ A >) ( 33 ) 



12 All these probabilities naturally appear if we postulate that the characteristic algebra is a nuclear one as we will do in the 
rigorous treatment of the subject elsewhere. 
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Then calling ( |§ eq. (V.1.7)) 

F { A %{z) = wp[Ba z {A)] 



G%{z)=wp[ az {A)B] (34) 



where z = t + 17, we obtain 



(I\ e -P H Be iHz Ae- iHz ) _ {I\Be lHt e~~< H ' A e - iHt e-^-^ H ) 
F a,b\z) - (/| e -/MT) - (I\e-f> H ) (35) 

But from eq., (pi| ) 

Tre _/Mr = (/|e -/,ir ) = / e-^dw (36) 
which is only convergent and for > 0. Analogously, the r.h.s. of eq. (|35|) reads 

J J B{oj,uj')e w,t e-^' ' A{J \tu)e^ UJt e^^-^ UJ du;duj' / J e'^du) 

Thus in order that F^ B {z) be analytic in z these integrals must be convergent and therefore the three real parts of 



the exponents: 0, 7 and, — 7 must be positive. Then the operator in (35) is analytic for 

< 7 < (37) 

Also from (||) 

G^ B (t) = F^ B (t + 0) (38) 

and both functions are analytic for the z = t + ij satisfying condition (|3^). Therefore KMS condition is satisfied. 

We have proved that the thermal equilibrium state exists in the system but we have not proved if the system 
spontaneously goes to this equilibrium (it "thermalizes" ) and we have not said under what conditions this phenomenon 
takes place. The reason is that the two powerful devises we have used: the diagonal form of hamiltonian H and 
information theory allow us to reach the above conclusion but they yield to a static and already thermalized diagonal 
of the density matrix. To study thermalization we must consider at least two subsystems and their interaction fce/orePI 
diagonalization and show how they reach the same temperature. In this case the method can be also used to show 
that, if the system is endowed with an adequate interaction, it naturally evolves to state (Eg), i- e. it "thermalizes". 
Thus the kind of interaction defines if the system thermalize or not. In fact, we have proved in paper || that the 
thermodynamic limit can by formalized with our method and we have developed a simple model with linear interaction 
(Fricdrichs model) where an oscillator is thermalized by a bath for small interaction. The case of big interactions will 
be treated elsewhere. 

Let us finally remark that a KMS state has a trivial explicit definition in a box. But in an unbounded space it was 
substituted by a state that satisfies KMS conditions with no explicit expression jl6| . Our method has allow us to find 
this explicit expression, namely 



III. THE CLASSICAL STATISTICAL LIMIT. 

In this section we will use the Wigner integrals that introduce an isomorphism between quantum observables O 
and states p and their classical analogues O w (q,p) and p w (q,p) Ej: 



3 Or after diagonalization but in an indirect way as in |l4| or jl5| 
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O w (q,p) =Jd\{q- ~\0\q + ^> exp(iAp) 

p w (g,p) = — L- / dA(p||g + A)(«-A|) exp(2zAp). (39) 



It is possible to prove that / dqdp p w (q,p) — (p\I) — 1, but in general p 14 ' is not always non negative. It is also 
possible to deduce that 



[p»\O w ) = J dqdpp w (q,p)O w (q,p) = (p\0), (40) 

and therefore: to the mean value in the classical space corresponds the mean value in the quantum space. Moreover, 
calling L the classical Liouville operator, and L the quantum Liouville-Von Neumann operator, we have 

L [p W (q,p)} = [hp] W (q,p) + 0(H/s), (41) 

where L p w (q,p) = i {H w (q, p) , p w (q, p)} pB , s is an action with the characteristic dimension of the system, and 

(hp\0) = (p\[H,0]). (42) 

Finally, if O = 0\0^ where 0\ and O2 are two quantum observables, we have 

O w (q, P ) = Of{q,p)Of{q,p) + 0{h/s). (43) 

We will prove that the distribution function p^(q,p), that corresponds to the state functional via the Wigner 
integral is a non negative function of the classical constants of the motion, in our case R obtained from the corre- 
sponding quantum operators H, Pi,..., Pn ■ 

From eq. ( |l4| ) we have: 



/>oo 

W lim p(t) = y2p r (uj )(u} 0l rr\ + V] / dwp r (u)(w,rr\ 
t ~*°° ^ Jo 



(44) 



so we must compute: 

P Z(q,p) = tt^- 1 J (w,rr||g + X)(q - X\)e 2ipX dX (45) 

We know from ||| section II. C, (or we can prove directly from eqs.([l4|-|l7])) that 

(u> ,rr\H n )=u>l (u, rr\H n ) = u n , 

(wo,rr|i*)=r?, (w,rr|i?) = r?, i = l,...,N (46) 

for n = 0, 1, 2, ... Using the relation ( pis] ) between quantum and classical products of observables and relation ( |40| ) 
between quantum and classical mean values, in the limit h — > (precisely when s — > 00, i. e. when the dimension of 
the system are very large compared with tif^\ we deduce that the characteristic property of the distribution p^ r (q,p), 
that corresponds to the state functional (u>, rr\, is^|: 



pZ{q,p)[H w {q,p)} n dqdp = u n , pZi<l,P)[Pi {q,p)] n dqdp = r? , (47) 



14 We can prove that H w (q, p), P™ (q,p),..., P™ (q,p) are constants of the motion using Heisenberg version of Liouville equation 
and eq. @. 

15 We will consider that we always take this limit when we refer to classical equations below. 

16 To simplify the demonstration we will consider that all the spectra are continuous as explained in footnote 7. 
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for any natural number n. Thus p^f r (q,p) must be the functional^ 

pZ(q,p) = S(H w (q,p) - Lo)5{P^{q,p) - r x )...5{P™ (q,p) - r N ). (48) 
For the distribution p^ or (q,p) corresponding to the state functional (ojo,rr\, we obtain 

pZAI'P) = 5{H W {q,p) - «o)5(P^(g,p) - n)...6(Pg(q,p) - r N ). (49) 
Therefore, going back to eq. (^) and since the Wigner relation is linear, we hav^j 

poo 

. Jo 



(50) 



po(H W (q,p),P^( q ,p),...,P^(q,p)) + p(H W ( q ,p),P^( q ,p),...,P^(q,p)) 

Also we obtain p*(q,p) > 0, because p r (wo) an d p r (u>) are non negative. 

Therefore, the classical state p^(q,p) is a linear combination of the generalized classical states pQ r (q,p) (where 
is either uj or lj), having well defined values f2, r\,..., of the classical observables H w (q,p), P x w (q,p),..., Pff(q,p). 
The corresponding classical canonically conjugated variables are completely undefined since neither pQ r (q,p) nor 
p^(q,p) is a function of these variables. So we reach, in the classical case, to the same conclusion than in the 
quantum case (see end of subsection II. A. 2). But now all the classical canonically conjugated variables do, Gti, ajv 
do exist since they can be found solving the corresponding Poisson brackets differential equations. We can also expand 
the densities given in eqs. ( ^8| - |5C| ) in terms of classical motions as shown in JsJ in great detail. In fact, as the momenta 



17 Let us precisely define the S of the next equation. In eq. (^2) we can make the canonical transformation q,p — > Q w , P v 
then it reads 



Then: 



where 



/ pH r {QT,pr)[pr\ n dQTdpr = r ?, 



N 

P Z{Q w ,p w ) = (v Q )- l i[6{py 

i=0 



V Q = / dQ w 



i = 0.1, ...N 



is the volume of configuration space that we will consider as bounded for simplicity. Nevertheless the case Vq — > 00 can be 
studied with the techniques of pap er Q . 
Then to prove eqs. ( [iq ) and ( |49| ) we can first write eq. as 



P* (<?,£>) = / dp(R)p R (uj )p uloR (q,p) + 



I MR) So 



where we have used eq. (|l5|). 
18 We will call: 



po(w,ri, rjv) = p r {oj ) = Pr 1 ,...,r JV (^o) 



p(uj,n,...,r N ) = pr(u) = p ri rjv (oj) 
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H w , P™ , P]q, or any function of these momenta, that we will call generically II, are also constant of the motion, 

then we have -^11 = —dH/da = 0, where a is the classically conjugated variable to II. So H is just a function of the 
II and: 

d dH(U) 

— a = = a7 ( 11 ) = const - ( 51 ) 

so: 

a i (t)=ro J (n)t + o j (0), j = 0,1,..., Ar. (52) 

Thus (going back to the old coordinates q,p) in the set of classical motions contained in the densities Jig ) and j4S| ) 
the momenta H w (q,p), P^ (q,p),..., Pj[(q,p), are completely defined and the origin of the corresponding motions, 
that we will respectively call ao(O), ai(0),...and ajv(0), are completely undefined (since /Q* does not contain 

these variables), in such a way that the motions represented in the last equation homogeneously fill the surface where 
H w , Pi and Pff , have constant values. If the system is bounded and integrable (see below), these surface turns 
out to be the a usual torus of phase space. This is the interpretation that we give to the density ( p0| ) which is just 
a function of the variables H , P± Pjy , but it is not a function of the classical conjugated variables do, ai,..., 
a at. The classical motions described by eqs. (^T|) and ( |j2"| ) will be used to characterize the dynamical system in the 
next section. Precisely, using these motions the evolution of the p w , close to equilibrium p^ can be interpreted as 
Frobenius-Perron operators |lj|[^} 
In conclusion. 

i. - We have shown that the quantum state functional p(t) evolves to a diagonal state /?*. 

ii. - This quantum state p* has its corresponding classical density p^(q,p). 

iii. - This classical density can be decomposed in classical densities that correspond to sets of classical motions 
where H w , Py Pff remain constant. These motions have H w (q,p), P™ (q,p),..., P^(q,p) — const and initial 
conditions ao(0), a± (0) , ajv(0) distributed in an homogeneous way. 

iv. - From eqs. (|4^-p0[) we will obtain thatj^J 

P 7(q,p)=p(H W (q,p),P^(q > p) ) ...,P^(q,p))>0 (53) 



IV. RELATION BETWEEN CLASSICAL STATISTIC, DYNAMICAL SYSTEMS AND 

THERMODYNAMICS . 

In this section we will characterize the classical motions we have obtained in the ergodic hierarchy and we will make 
contact with classical thermodynamics. 



A. The microcanonical ensemble. 



1. The notion of isolating and non-isolating constant of the motion. 



The classical constants of the motion H w (q,p), P^ v (q,p), Pff(q,p) can be rigorously classified as |2(J, VM , P2| : 
L- Global or isolating constant of the motion, that we will call " H" , precisely: 

H™(q,p) = H w (q,p), Hf(q,p)=P^(q,p),...,Hf(q,p)=PY(q,p) (54) 

when the conditions 

Hj v = r i , i = 0,..., A (55) 



19 This is an important step towards classicality since the system has lost its quantum characteristic (e. g. totality, con- 
textuallity and non-locality according to |l9|) and we can resolve the classical density matrix evolution in a set of particle 
motions. 

20 From now on we will forget the bound eigenvalue u>o- 
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(where r$ — f2) define global sub manifolds (tori in the bounded case) M(ro, ■■■, r A)i of phase space where the 
trajectories necessarily move, for each set of constants (ro, The dimension of M. is 2(N + 1) — (A + 1) = 2N 

-A + l 

ii.- Local or non isolating constant of the motion, that we will call: 

J^(q,p) = P^ +1 (q,p), J^ A (q,p) = P%(q,p) (56) 

when the conditions 

J? = r j+A j = l,..,N-A (57) 

do not define any global sub-manifold, since they are just local The jj^ can only be considered as local coordinates 
on the manifold M(ro, ta)- 

When A = N we say that the system is integrable, when A < M we say that the system is not integrable. Let us 
consider both cases. 



2. Integrable systems. 

In this case all the P are isolating constants of the motion (H) and there are no J. Then the situation is like the 
one described at the end of section III: condition ( |55| ) foliates phase space with submanifolds A4(ro, rjy), labelled 
by the constants (ro, ...,rjv). This submanifold would be tori if the system is bounded. As we have already said we 
will only consider the bounded case^j. In the tori the motion of the configuration variables is given by eq. j5^). 
In the generic case the Wj(H) are not rationally dependant (or non-commensurable). Then the trajectories of the 
configuration variable fill each torus in a dense way. Therefore the motion is ergodic in each torus. Moreover, we can 
independently see from eq. ( |53| ) that there is a unique equilibrium state in each torus 

pY(q,p) = p(ro,n, -,rjv) (58) 

which is constant in the submanifold. Therefore we have a microcanonical equilibrium in each torus. 

At this point we know that there is a unique stationary equilibrium state p^(q,p) = p{tq, rjy) but we do not 
even know if the evolution converges to pjf . So now we will prove, using the quantum equations ( DTI) or (44), that 
when t — > oo, any p w — ► in a weak way. 

In the quantum case we have 



lim (p(t)\0) = (p*\0) (59) 



Now, quantum products can be changed into classical products according to eq. (BO) so: 



and we have proved that: 



lim(p w (t)\O w ) = (pW\O w ) (60) 

t— *oo 



W lim p w (t) =pf (61) 

t — >oo 



Therefore any p — » pj; weakly when t — > oo. 



21 Isolating constant of the motion would be the "simple" constant of the motion in [^| p. 60. 
22 Let us list the introduced dimension: 

i. -The total dimension is 2(N + 1). 

ii. -The number of the isolating constants is A + 1. 

in. -The number of the non isolating constant si TV — A. 

iv. -The number of configuration coordinates is N + 1. 

v. - The dimension of M is 2/V - A + 1 

23 The unbounded case requires a new singular structure in the states (see 0). This generalization will be studied elsewhere. 
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3. Non integrable case. 



In this case not all the P are isolating constants of the motion so there are H and J. Then the trajectories must 
be dense in a domain 2Xjo, ■■■,va)c M(r , ■ ■■,ta) of dimensions 2N — A + 1. If not the dimensions of T>(r , ...,r^) 
would be < 2N — A + lrj. In this case a new global constant must exist and there would be A + 2 global constants. 
But this is impossible since A + 1 is the total number of these constants. It is quite clear that, as the equilibrium 
classical density must be globally defined in T>(ro, ta) the (q,p), ■ J]^_a{<1,p) cannot be explicit variables of 
P^fefOEB so we mus t just have: 

p™(q,p) = p(H™(q,p), HF(q,p),...,H%(q,p)) (62) 
But on 2?(ro, Ta) we have that: 

H w (q,p) = ro = const., H^ v (q,p) = r± = const., H^ (q,p) = Ta = const. (63) 
so on V(tq, ta) it is: 

PV(Q,P) = P( r o, n, ...,r A ). (64) 

and we find the unique equilibrium for each T>{r§, ...,ta)i namely for each set of constants (r*o, ...,ta)- Now we can 
follow the reasoning of the previous section. Phase space is now foliated by submanifolds M(tq, ta) of dimension 
2N — A + 1. The changes are that now not all the coordinates of these submanifolds are configuration variables^, 
there are also momentum variables, therefore we cannot use the reasoning about the zuj(ll) not rationally related in 
this case. Nevertheless in each T>(tq, ...,Ta)C M(ro, ...,Ta) there is a unique equilibrium state ( p4| ) so using theorem 
4.3 of ref. |Q we conclude that the motion is ergodic in 2?(ro, ...,7^4). Now we can repeat the reasoning of eqs. ( |59| ) 
to ( |6l| ) to show that the equilibrium (Q) is weakly reached. Since this equilibrium is a constant we again find a 
microcanonical ensemble. 

But, at this point we may ask ourselves why function p(H^ (q,p), H™ (q,p), (q,p)J™ (q,p), jff_ A (q,p)) 
looses its J variables. To explain this fact we may say that really space O must contain physical measurable ob- 
servables. But only the CSCO can be measured in an independent way since their observables commute. Moreover 
as the classical momenta (q,p), jff_ A (q,p) have an ergodic motion so it is reasonable to consider that the 
quantum analogues J±, Jn-a, cannot really be measured, even at the quantum levej^. Then the set {Hq, Ha} 
is the relevant measurable CSCO, and it is only possible to measure the H,H\, ...,Ha- Then the isolating constant 
{Ho, ...,Ha} remains as the only characters in the quantum or classical play. 

Phrased in another words the J variables can only be diagonalized locally, while the diagonalization procedure of 
section II. A. 2 was thought as a global one. So it is better to consider that the unitary operators U of eq. (|l^) only 
diagonalize the indices of the H (that we will call r) and do not diagonalize the indices of the J (that we will call m). 
Then we obtain a basis where the coordinates of the stases read ( see ([l3|)) 

p{io]rmr'rn' — Prmm' '\W)$rr r (^^) 

Then eqs. (Elf) and (jig) become 



24 We may say that ^4 + 1 tori are not broken and that A tori are broken. 

25 From its definition, in eq. ( |l3| ) or in eq. (|5C)|), p is just an ordinary function of global variables. Therefore it cannot be 
considered as a function defined using local coordinates. 

26 Regarding the configuration variables it is clear that, in the non-integrable case, none of then is a global constant of the 
motion, that further reduces the dimension of M (or T>). In fact: 

L- The preserve tori satisfy an irrationality condition ( ^l), eq. (3.4.12)) so the corresponding ratios of the frequencies are 
irrational and the trajectories are dense in these tori. 

ii.- In the broken tori the trajectories are chaotic. 

Nevertheless, if in a particular case there is a configuration variable X that turns out to be a global constant of the motion it 
can be considered among the "if". Then we will work essentially in the manifold X = const, and nothing will change. 

27 We can only measure dynamical variables when they can be considered as constants in time, at least in the period of 
measurement. If a constant is not global it is only constant in time in a local coordinate system, i. e. it is not really physically 
constant. 
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p* = W lim p{€) = } \ dujp rmm i{ijj){uj,rmrm'\ (66) 

t^oo J 
rram' 

and 

Pi = E J du J P t rmm ,\urm)(uorm l \, i = l,2,...,A (67) 

rram' 

and so on for the rest of the equations. 

In order to follow this analysis we must "trace away" the m since the corresponding observables cannot we measured, 
i. e. they cannot be considered classical in a global way. Essentially we must consider the J operators and the m 
indices as inexistent in space O. As O is the space of all measurable observables this fact must be considered as the 
following change in the observables of eq. (||) 

O(oj) (68) 

(we consider only the diagonal term since we are only concerned in the classical part and we neglect the luq term 
since we are systemically forgetting the ground state) . In this way the m index has a " spherical symmetry" and they 
measure nothing. Then the relevant part of eq. (||) reads 

(O) p = J dujp(u}) rmr ' m >0(uj) rmr ' m > — J dujp(uj) rmr ' m '0(uj) rr iS mm ' = 



rmr rn rmr m 



rr' \ rn 



u)rmr'mj O(to) rr > (69) 

Calling Y^tm P(v)rmr'm = p(<*>)rr' we obtain the "traced" equation: 

{0) p = J duip(uj) rr i O(oj) rr i (70) 

rr' 

and the m indices have disappeared. From now on we can work with only the r indices (ro, ta) and in this way 
the p of eq. (^8|) becomes the p of eq. (|64|) solving the problem. 
Two comments are in order: 

L- From what we have said eq. ( [l4| ) can be considered as the quantum version of the microcanonical equilibrium. 
Therefore the conditions to obtain decoherence are equal to those to obtain microcanonical equilibrium states. But 
this state will be different according to the number of J. 

ii.- H has a continuous spectrum, so this is also the case of L. (see (fl^). Thus, albeit the 0(h/s) of eq. (fH]), we can 
conclude that the classical L has a continuous spectrum. Moreover, normally quantum corrections make discrete the 
classical continuous spectra, and not viceversa. Even more, the spectrum of classical Liouville operators are usually 
continuous. Then L, most likely, has a continuous spectrum. This is one of the characteristic properties of mixing 
(and therefore ergodic) systems, the spectrum of their evolution operators are continuous. This is another way to see 
that the flow is mixing in T>(ro, ...,r^) p3{, so the existence of a weak limit is natural. 



B. The canonical and grand canonical ensemble. 

We have proved that our system reach a microcanonical equilibrium. Then, the frequent presence of canonical 
equilibrium can be explained by at least three reasonings: 

i. - The thermodynamic limit. In this limit microcanonical and canonical densities coincide (see e. g. p4|). 

ii. - Canonical subsystem of a microcanonical system. A small subsystem of a big microcanonical system 
(that can be considered as a thermic bath) is canonical p4j]. So, if the only isolating constant of the motion is the 
energy H, we have reach to the notion of canonical equilibrium: 

p w „ e -w w (71) 
where (3 = T _1 . In the more general case of section III we would arrive to the conclusion that (cf. p5|) 
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p? ~ eM-PK ~ 7i^r - - - 7A#I ) (72) 

namely generalized grand-canonical ensembles. 

iii.- Information theory. If we only use unbiased distributions as in section II. B the exponential will naturally 
appear |l3j. Moreover our CSCO are numerable set, even if some of its observables have continuous spectrum, so we 
can use the formalism of the quoted paper []l3|| , based on the information theory, with only one modification: to use, 
for the observables with continuous spectrum, the classical continuous version of Shannon (eq. (pl|)). Then we 

will again find eqs. ( |7l| ) or (|72]). 

In all these ways we can obtain the canonical ensemble, to define temperature and to begin the development of 
classical thermodynamics. As at the end of section II. B. 2 we can say that the system will thermalize or not according 
to the kind of its interaction. 

We have a final verification: to relate the equations of II. B with the classical equations above. 

Let us begin with the "trace of an operator" eq. (^4j) in the simplified CSCO {H }. The classical object corresponding 
to (uj\ is by eq. ( |4§| ) 

p™(q.p) = 6(H w (q,p)-u>) (73) 
so the classical object corresponding to (I\ is 

I W (q,p) = J 6(H w (q,p) - u)du = 1 (74) 

Therefore the classical trace of an operator corresponding to (I\0) is 

(I W \0 W ) = J W (q,p)dqdp (75) 

namely the integral usually associated with the trace. Finally the classical equation corresponding to (Eq) reads 



W ^ A >- fe-^ w Mdqdp {7b) 

i.e. eq. (|7ll), if we consider pY as a functional of A and for the simplified CSCO {H}. In the general case the 
functional version of ([72]) would be 

w f as Je-^ W -^ V --^A^(q,p)dgdp 

w e w - j-m^w-^^iiojqTp [ ] 



i. e., the classical version of (E9h restricted to the isolating members of the CSCO. 



V. LOCALIZATION. 



The fate of a classical statistical system is to remain unlocalized in the phase space or to localize around a classical 
trajectory, in which case the system becomes classical. This fate depends in the potential acting in the system 
(which can be localizing potentials or not) and the initial conditions. Let us consider a classical distribution p(q,p) 
with support of volume AV^ N+1 \ It would localize if AV^- N+1 ^ — > when t — > oo. But this is impossible since we 
are dealing with an hamiltonian system where, according to Poincarc theorem, the volume of phase space remains 
constant. Nevertheless we see localized classical objects. This is only possible of the system corresponds to those 
studied in section IV. A. 3 where there is a volume AVq A+1 \ that corresponds to the projection of the support 
of p{q,p) over the subspace of the observed dynamical variables H, that vanishes when t — > oo, while the volume 
AVy N A ^ corresponding to the unobservable J dynamical variables diverges when t — > oo being 

AV<j N - A) .AV^ A+1) ~ AV( N+ V = const. (78) 

Therefore the subsystem of the observable variables is not hamiltonian and the Poincare theorem is by-passed. We 
conclude that only non integrable system become classical and that the "non-observability" invoked in section IV. A. 3 
is essential to obtain a final classical mechanical limitp| 



28 We will give examples of the localization phenomenon elsewhere. 
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VI. CONCLUSION. 



In usual text books of thermodynamic or statistical mechanics the ergodicity of the system or the micro-canonical 
state are postulated. With our method we have proved these postulates for systems with evolution operator endowed 
with continuous spectrum onlyp^|. Is this a limitation of our result? Quite on the contrary, it characterizes the systems 
where thermodynamic and statistical mechanics can be used: ergodic or mixing systems, being the latter one endowed 
with continuous evolution spectra p3f . Of course if we work on a box the evolution operator may have a discrete 
spectrum and we cannot reach to this conclusion. But we know that in this case sooner or later we must expand 
the box to infinity. Then, either we use hand waving argument: e. g. that the distance between the values of the 
spectrum is very small etc., or we use the rigorous result that in the limit the spectrum is continuous, and from this 
fact deduce we the physical properties of the system, as we have done in this paper 
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